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We propose a method to implement cavity QED and quan- 
tum information processing in high-Q cavities with a single 
trapped but non-localized atom. The system is beyond the 
Lamb-Dick limit due to the atomic thermal motion. Our 
method is based on adiabatic passages, which make the rel- 
evant dynamics insensitive to the randomness of the atom 
position with an appropriate interaction configuration. The 
validity of this method is demonstrated from both approxi- 
mate analytical calculations and exact numerical simulations. 
We also discuss various applications of this method based on 
the current experimental technology. 

PACS numbers: 03.67.-a, 42.50.-p, 42.50.Gy 



I. INTRODUCTION 

Trapping of single atoms in high-Q cavities opens up 
exciting possibilities for observation and manipulation of 
the dynamics of single particles and for control of their 
interactions with sing le-mode photons f|[||j|. Such 
possibilities could have wide applications, such as for gen- 
eration of non-classical or entangled optical pulses ||[7j , 
for observing strong cavity- QED effects j||5|||, and more 
remarkably, for implementation of quantum communica- 
tion and computation [j9|-fl3[ . The trapping potential for 
confining single atoms can be created by diverse avenues, 
including by the cavity QED light itself [H||, by an ad- 
ditional far-of-resonant trapping (FORT) beams a, and 
by combining single trapped ions with high-finesse op- 
tical cavities Jl4],[l5]]. In this paper, we will direct our 
attention principally to trapping in cavity QED by way 
of an additional FORT beam, although our results are 
applicable to broader settings. 

The first experiment to achieve strong coupling in cav- 
ity QED with trapped atoms was that of Ref . [g] , which 
employed an intracavity FORT beam and reported trap- 
ping lifetimes of 28ms. By now, this experiment has 
attained much longer trapping times, with recent work 
demonstrating lifetimes in excess of Is |||n|. By con- 
trast, atomic localization by way of the cavity QED field 
itself has led to a trapping within a single axial well with 
mean trapping time r 340/is Q and to localization 
across many axial wells with mean time r w 280/is @j . 

The long trapping times achieved with an intracav- 
ity FORT beam set the stage for diverse applications in 
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Quantum Information Science, which motivates the cur- 
rent analysis. However, one of the main obstacles to the 
experimental demonstration of these applications is that 
the position of the trapped atom is not well fixed within 
the cavity. The coupling rate g between the atomic in- 
ternal levels and the cavity mode depends on the atom 
position r through the relation 

9 (r) = goX (r) (1) 

with the mode function 

X (r) = sin (k z) cxp [- (x 2 + y 2 ) /w 2 ] , (2) 

where go is the peak coupling rate, wq and k$ = 2tt/Xq 
are respectively the width and the wave vector of the 
Gaussian cavity mode, and z is assumed to be along the 
axis of the cavity. Due to the randomness of the atom 
position r, we have an unknown randomly changing cou- 
pling rate g (r) . Most of the applications of this setup 
assumed a fixed known coupling rate g. Therefore, be- 
fore the experimental demonstration of these schemes, 
first one needs to solve the problem associated with the 
random coupling. 

Intense experimental efforts have been taken to local- 
ize the atom inside the cavity so as to fix the coupling 
rate g(r), with notable recent success attained via ion 
traps p^ , p^ j. In the cavity QED experiments employing 
cold atoms and without FORT beams |l|, [i~7j , p^ | , atoms 
were dropped through the cavity and followed random 
trajectories with large axial heating. As a result, the 
magnitude and the sign of g (r) were not well controlled. 
With a FORT beam and with current experimental ca- 
pabilities HQ, an atom can be trapped inside one 
potential well along the cavity axis with a fixed sign of 
g (r) . But the atom still has appreciable kinetic energy 
and is not fully localized, leading to significant variations 
in the magnitude of the coupling rate g (r) . 

The randomness of the coupling rate g (r) comes from 
several contributions: first, the trapped atom is still quite 
hot in the current experimental setup. Its kinetic en- 
ergy from the thermal motion is typically lower but not 
much lower than the depth of the trapping potential. The 
atom's oscillation amplitude d in the trap is comparable 
to the optical wave length Ao, so it does not satisfy the 
usually assumed Lamb-Dick condition d <C Xq. Due to 
the thermal motion of the atom, the coupling rate g (r) 
typically has a variation within a fact of 2 with the cur- 
rent experimental technique. Certainly, the atom will be- 
come better localized as cooling techniques are adapted 
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to cavity QED and its energy is reduced | J19|J2C| ] . How- 
ever, due to the presence of the cavity and the trapping 
potential, it is still experimentally hard to achieve ef- 
ficient cooling inside the cavity fll9|j20|| . Furthermore, 
even if we assume that the atom has been pre-cooled 
and localized initially to the Lamb-Dick limit, the imple- 
mented application protocols will still tend to heat the 
atom due to photon recoils from the spontaneous emis- 
sions pi 22 1 . As a result of the heating, the atom may go 
out of the Lamb-Dick limit after a short time. Finally, 
even if we neglect all the motional and the heating ef- 
fects of the trapped atom, there is still some uncertainty 
of the coupling rate. The intracavity field of the FORT 
beam forms many potential wells inside the cavity, and 
in current experiments, one can not control and does not 
know precisely in which well the atom is trapped. The 
FORT beam has a wavelength Xp different from the cav- 
ity QED wavelength Ao, so, even if the atom is kept to be 
very cold and well localized at the bottom of the trapping 
potential well, we still might not know exactly the cou- 
pling rate since the bottoms of different potential wells 
have different coupling rates |22j ]. 

Here, to overcome these difficulties, we propose a 
method to do cavity QED and quantum information pro- 
cessing directly with hot atoms with an inhomogeneous 
distribution in position and/or time varying locations. 
The method is based on adiabatic passages with an ap- 
propriate interaction configuration. Normally, schemes 
based on adiabatic passages are more insensitive to cer- 
tain parameter changes compared with the corresponding 
Raman schemes. Some initial indication of insensitiv- 
ity of the adiabatic passage scheme to certain parameter 
changes was already illustrated in |2^| for some cavity 
QED scheme. However, to make the whole system dy- 
namics insensitive to variations of the coupling rate g (r), 
we also need to design an appropriate interaction configu- 
ration. The relevant dynamics of adiabatic passages are 
determined by the relative ratio between different cou- 
pling rates, and are almost independent of their absolute 
values. Thanks to this property, with an appropriate 
design of the interaction configuration, we can make dif- 
ferent coupling rates have the same dependence on the 
atom position r, and therefore the system dynamics de- 
termined by their relative ratios will become independent 
of r. As a result, though the atom position may be un- 
known and time dependent, the output signal from the 
cavity is still well controllable and has definitely known 
properties. Note that the method described here is differ- 
ent to some previous quantum computation schemes with 
hot trapped ions J2^j2^] , where the Lamb-Dick condition 
is still required. 

The paper is arranged as follows: In Sec. II, we ex- 
plain the basic idea of the method, and then describe 
and solve the model Hamiltonian analytically within the 
adiabatic approximation. In Sec. Ill, we give an exact 
numerical simulation of the model, with the emphasis 
on checking the validity of the introduced approxima- 
tions and calculating various kinds of noise magnitudes 



relevant for the on-going experimental efforts. The cal- 
culations show that we can get reasonably good signal- 
to-noise ratios with typical experimental values for the 
parameters. In Sec. IV, we discuss various applications 
of this method. By incorporating this method into some 
previously-known schemes, we show that with hot non- 
localized atoms, one can still realize many kinds of cavity 
QED and quantum information processing schemes, in- 
cluding, for instance, the controllable single-photon or 
entangled photon source, quantum communication be- 
tween cavities, atomic entanglement generation, telepor- 
tation, and Bell inequality detection. Section 5 gives a 
synopsis of parameters relevant to our current experi- 
ment for a single atom trapping with a FORT beam at 
Caltech P,^|Jlq|. We summarize the results in the final 
section. 



II. CAVITY QED WITH A NON-LOCALIZED 
TRAPPED ATOM: THE SCHEME 

A. Basic idea 

First, we explain the basic idea of this method by con- 
sidering a single trapped atom, which has three effective 
levels \g), |e), |s), as shown in Fig. 1. The two ground 
states \g) and |s) can correspond, for instance, to sub- 
Zeeman levels in the F = 3 and F = 4 manifold respec- 
tively for the cesium atom. The transition |e) — > \s) is 
coupled resonantly to the cavity QED mode a with a cou- 
pling rate g (r) in the form of Eq. (1). A classical laser 
field e (t) incident from one mirror of the cavity (see Fig. 
1) drives the transition \g) — ► |e) through another cav- 
ity mode a'. We assume for simplicity that a and a' have 
the same spatial mode structure with the same frequency 
(for example, they can be of different polarizations). The 
driving laser s (t) is resonant to the transition \g) — > |e), 
so it is far-of-resonant to the cavity mode a' with a large 
detuning uj gs , where uj gs denotes the splitting between 
the levels \g) and |s). Due to the off-resonant driving by 
s(t), a' can be described classically by its mean value 
(a') = a (t) e~ tuJ !><! t ( U!ge [ s the frequency splitting be- 
tween the levels \g) and |e)), which couples resonantly to 
the transition \g) — > |e) with a Rabi oscillation frequency 
Q (r, i). Since a and a' have the same spatial mode struc- 
ture, the Rabi frequency fi (r, t) will depend on the the 
atom position r by the same mode function x( r )i i- e -i 
f2 (r, t) can be factorized as 



ft (r, t) = n (t) X (r) = r o g a (t) X (r) 



(3) 



where r Q represents the fixed ratio of the Clebsch-Gordan 
coefficients for the transitions \g) — > |e) and \s) — » |e). 
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FIG. 1. Schematic setup. Left side: a single atom trapped 
in a high-Q cavity, which is driven by a classical laser pulse 
e (t). Right side: the relevant atomic level structure. 

To understand the basic idea of this method, let us 
first look at a very simplified picture by neglecting the 
coupling of the mode a to the cavity output. The system 
is then described by the following simple Hamiltonian in 
the rotating frame (setting % = 1) 



H sim = SI (r, t) a eg + g (r) aa es + H.c. 



(4) 



where = (v\ (/i, v = g, e, s) are the atomic tran- 
sition operators, and H.c. stands for the Hermitian con- 
jugate. The Hamiltonian i? s im has the well-known dark 
state \D) (the instantaneous eigenstate with a zero eigen- 
value) with the form p5l 



\D) = 



1 



ls(r)| 2 + |fi(r,i)| 
1 



GKr)| 5 )|0>-n(r,t)|*>|l>) 



(| S )|0)-r o a(t)|a)|l», (5) 



where |0) and |1) represent the zero and the one photon 
state of the cavity mode a. Note that the dark state \D) 
actually only depends on the ratio between the parame- 
ters g (r) and fi (r, t) , so it becomes independent of the 
random atom position r with the interaction configura- 
tion specified above. If we start with the atom in the 
ground state \g), and gradually increase the Rabi fre- 
quency fl (r, t) , under the adiabatic approximation, the 
state of the system will remain in the dark state 
which gradually evolves into the final state \s) |1). Due 
to the independence of the state \D) on the variable r, 
the relevant dynamics of this adiabatic evolution also be- 
comes independent of the random atom site r. This is 
the basic idea of the method to eliminate the influence of 
the randomness in the coupling coefficient g (r) . 

Note that to make the dark state and the relevant dy- 
namics independent of the random atom position r, the 
driving pulse and the cavity mode need to have the same 
spatial mode structure. This is why the classical driv- 
ing pulse is matched to spatial mode of the cavity field, 
both along the cavity axis and transversely, which is rou- 
tinely accomplished by way of illumination from one side 
mirror of the cavity. This configuration is different from 
the original proposals for adiabatic dynamics in cavity 
QED p5j in which the propagation direction of the driv- 
ing pulse is perpendicular to the cavity axis with uniform 



illumination intensity. It is also distinct from the config- 
uration employed in some recent interesting experiments 
directed toward achieving a single-photon source 0,|l8| , 
which likewise employed uniform illumination transverse 
to the cavity axis and for which the atom is not localized 
axially. As a result, in these experiments some of the dy- 
namics, such as the output pulse shape and phase, still 
depend on the unknown position of the randomly moving 
atom, and are thus not fully controllable and reversible, 
as has been seen from the experiments. 

To guarantee an adiabatic evolution, we need to ful- 
fill the adiabatic condition, which means the evolution 
time T should be significantly longer than the frequency 
gap 5 between the dark state and some other eigenstates 
of the Hamiltonian H slm . The error probability due to 
the non-adiabaticity is estimated by p a d = 1/ (ST) 2 . For 
the Hamiltonian H s i m , the frequency gap S is given by 



5 = y \g (r)\ + |fi(r, t)\ . Thus, the adiabatic condition 
\g(r)\ 2 + \Q(r,t)\ 2 T 2 > 1 depends on the atom po- 
sition r. If the coupling coefficient g (r) changes by a 
factor of 2, the error probability p ac i will change by a fac- 
tor of 4 for the same evolution time T. However, if T is 
sufficiently long, the error probability p a d remains to be 
small, and the relevant system dynamics will be still very 
insensitive to the randomness of the atom position. To 
estimate p a d, we can use the average value of the coupling 
rate g (r). 

In the above simple picture, we neglect the coupling 
of the mode a to the cavity output. This is only a valid 
picture in the good-cavity limit with the evolution time 
T <C where k is the cavity decay rate. However, in 
practice, it is better to operate the system in the limit 
with T > 1/k. There are several advantages of operating 
the system in this limit: first, without the requirement 
T« 1/k, it is easier to satisfy the adiabatic condition for 
which T should be sufficiently long; second, in this limit 
it is easier to modulate the Rabi frequency fi (r, t) by 
changing the intensity of the driving laser e (t) incident 
from one side mirror of the cavity. In this way, one can 
efficiently control the pulse shape of the cavity output by 
modulating the shape e (t) of the driving laser, which is 
useful for many applications. In the limit T > 1/k, we 
need to take into account from the beginning the coupling 
of the mode a to the continuum cavity output, and the 
whole system will then have infinite levels. We will de- 
scribe in the next subsection this more involved interac- 
tion configuration. The above simple three-level picture, 
though it does not describe the real experimental config- 
uration, does help for the understanding of the basic idea 
of the adiabatic method. 
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B. Theoretical model and its approximate analytical 
solution 

Now we look at the more complicate theoretical model, 
which includes the coupling of the mode a to the contin- 
uum cavity output. If we adiabatically apply a classical 
diving pulse e (t) as shown in Fig. 1 , one photon will 
be emitted from the transition |e) — > |s), and the cav- 
ity will output a single-photon pulse. We want to show 
in the below that this single-photon pulse has a definite 
pulse shape which is independent of the randomness in 
the atom position r and in the coupling rate g (r) . In this 
way, although the atom position and the absolute value of 
the light-atom coupling rate are not fully controlled, we 
can nevertheless fully control the properties of the output 
single-photon pulse by modulating the driving laser pulse 
e (t). This is an important feature for many applications 
of this setup, which we will discuss in Sec. IV. There are 
several equivalent ways to describe the coupling of the 
mode a to the continuum cavity output. Since we want 
to calculate the output pulse shape within the adiabatic 
approximation, it is convenient to use the Hamiltonian 
approach |^8|,^9| . The whole Hamiltonian, including the 
coupling to the cavity output, has the following form in 
the rotating frame |2^] 

H = (A - ij s /2)a ee + [ft (r, t) a eg + g (r) aa es + H.c] 

dw [a)b{uj) - a6 t (ui)] (6) 



dui [ub' (to) b (w)] 



where b(uj), with the standard commutation relation 
[b (uj) , b^ (a/)] = 8 (uj — uj'), denote the one-dimensional 
free-space modes which couple to the cavity mode a. We 
only need to consider the free-space modes within a finite 
bandwidth [ui se — ojf,, ui se + uib] with the carrier frequency 
w se (uj se is the frequency splitting between the levels \s) 
and |e)), since all the modes outside of this bandwidth 
have negligible contributions to the dynamics due to the 
large detuning (larger than ujb). Within this bandwidth, 
the coupling between b (lj) and the cavity mode a is ap- 
proximately a constant, and we denote it by \J k/2tt for 
convenience, where n is the effective cavity decay rate, 
as we will see. The bandwidth u>b should be chosen to be 
much larger than k, but still much smaller than u> se . 

We have assumed that the driving laser and the cav- 
ity mode a couple resonantly to the corresponding free- 
space atomic transitions. However, we emphasize that 
our scheme still works for the case of off-resonant cou- 
pling. By considering the off-resonant scheme, there is 
no win with respect to losses due to the atomic decay, 
since in this case the time scale also slows down. So it 
suffices here to consider the resonant coupling case. How- 
ever, in the Hamiltonian (6), it is still helpful to include 
a single-photon-transition detuning A to account for the 



trapping potential difference for the levels \g) and |e) in- 
duced by the FORT beam (this potential is basically the 
same for the levels \s) and \g) for a FORT beam with 
linear polarization as in our current experiments). The 
potential difference between the level \g) and |e) in gen- 
eral depends as well on the random atom position r. 

The imaginary part of the Hamiltonian (6) accounts 
for the spontaneous emission loss, where 7 S denotes the 
total spontaneous emission rate of the upper level |e). 
In writing this form, we have assumed that the sponta- 
neous emission photon escapes and that the atom after 
a spontaneous emission will not be repumped. This is 
a good assumption for the interesting region where the 
spontaneous emission loss is not big, and the atom thus 
has a very small probability to be repumped after emit- 
ting a spontaneous emission photon. As a result of this 
assumption, the spontaneous emission only contributes 
to the leakage error which is properly represented by Eq. 
(6) @. 

We treat the atom position r in the Hamiltonian (6) as 
a classical stochastic variable, and neglect its quantum 
nature. This is a good approximation for the current 
experimental situation where the atom is still quite hot. 
There have been some analysis of the noise from quantum 
motional effects in high-Q cavities with very cold atoms 

& 

We start with the atom in the ground state \g), and 
then apply a classical driving pulse e (t) . This pulse 
can efficiently control the time evolution of the Rabi fre- 
quency VL (r, t) in the Hamiltonian (6). To see this, we 
write the input-output equation for the cavity mode al 



ai= —iujepdi - 



-al - 



(*) 



(J) 



where a[ n (t) is the field operator for the input driving 
pulse coupling to the mode a/, with (a[ n (t)) — e (t) and 

a[ n (t) , a^ n (t') = 6 (t — £'). By assumption, the mode 

al has the same frequency as the mode a which is res- 
onant to the free-space atomic transition \s) — > |e), so 
the eigen- frequency of al is uj se . Such a situation cor- 
responds, for example, to the case of the (a, al) modes 
of orthogonal polarization, but degenerate in frequency, 
although this is not an essential requirement. In Eq. (7), 
we have neglected the small depletion of a/ caused by the 
coupling to the atomic transition a eg since al is driven 
by a strong classical pulse e (t) which dominates its time 
evolution. We write the mean values of a' and a' in (t) as 
(a') = a (t) e - iw 9«* and (a[ n {t)) = e (t) = e(t) e^'* 1 , 
where s (t) is the slowly- varying amplitude of the driving 
laser. Form Eq. (7), we get a time evolution equation for 
the mean value a (t), which has the following immediate 
solution 

a(t)= [ e(T)e (4 ^- K/2)( *- T W. (8) 
Jo 

The variation rate of e(r) is characterized by the in- 
verse of the operation time T (the pulse duration) , which 
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is typically much smaller than the hyperfine frequency 
splitting Ljg S (about 9GHz for cesium atoms). Hence, a 
partial integration of Eq. (8) yields 



a(t) 



e{t) 



_ (iu) gs -n/2)t 



6(0) 



—iujgs + n/2 



l + o 



1 



UJgsT 



(9) 



We assume that s(t) gradually increases from zero with 
e (0) ~ 0. Then, within a good approximation, we have 
a (t) oc e(t) from Eq. (9). In the following, without loss 
of generality, we assume a (t) to be real by choosing an 
appropriate constant phase of e (t). The time behavior of 
the Rabi frequency (r, t) is completely determined by 
a (t) (note that (r, t) = r o g a (t) x (r) from Eq. (3)), 
that is, by the amplitude e(i) of the driving laser. 

The dark state (5) can be rewritten as \D) = 

cos%) |0) - sin0|s) |1), with cos0 = l/\Jl + \r a(t)\ 2 
independent of the atom position r. The state \B) com- 
plementary to the dark state is usually called the bright 
state with \B) = sm9\g) |0) + cos6» \s) |1). To solve the 
dynamics governed by the Hamiltonian (6), we can ex- 
pand the state | \&) of the whole system into the following 
superposition 

I*) = (c d \D) + c b \B) + c e \e) |0» ® |vac) + \s) |0) ® \(fi) , 

(10) 

where |vac) denotes the vacuum state of the free-space 
modes b (u>), and 



m) = 



dujcja' [uS) | vac) 



(11) 



represents the state (not-normalized) of the single-photon 
output pulse. The coefficients q, q,, c e and c u in Eq. 
(10) are time dependent. At the time t — 0, we have 
Cd = 1, Cb = c e = c u = and cos 8 = 1. After applying 
a classical driving pulse e(t), cos 9 slowly changes with 
a (t) , and we need to compute the time evolution of all 
the coefficients Cd, q,, c e , c w in Eq. (10) by substituting 
|\&) into the Schroedinger equation idt |\&) = H \^>). 

To go on with this task, let us first take the adia- 
batic approximation, which assumes the time derivative 
dt cosd w 0. As a result, dt \D) and dt \B) become negli- 
gible. We will check the validity of the adiabatic approx- 
imation and calculate various non-adiabatic corrections 
in the next section through numerical methods. In the 
adiabatic limit, the populations in the bright state \B) 
and in the excited state |e) are negligible, so we assume 
q, sa c e w 0. The coefficients Cd and c u satisfy the fol- 
lowing evolution equations 



- \/k/2tv sin 9 



jdcu, 



\/ k/2itc, 



d sm( 



(12) 



(13) 



Equation (13) has the solution 

CuJ (t) = f e^-^Cd (t) sinfl (t) dr, (14) 

Jo 

which, substituted into Eq. (12), leads to 
rt sin [w b (t - t)] 



Cd = - — sin 

Z7T 







t - T 



(k/2) c d sin 2 



Cd (t) sin (r) <ir 

(15) 



The approximation in Eq. (15) is valid since the band- 
width ujb satisfies lu^T 3> 1, where the operation time T 
characterizes the time scale for a significant change of Cd 
and sin 9. Therefore, the dark-state coefficient Cd satis- 
fies the cavity free-decay equation, with the decay rate k 
replaced by the effective rate K sin 2 9. This can be easily 
understood since sin 2 9 is the probability of the compo- 
nent |s) |1) in the dark state \D), and it is exactly this 
component that couples to the cavity output. Equation 
(15) has the straightforward solution 



Cd 



cxp J sin 2 9(r)dT 



(16) 



We want to know the single-photon pulse shape / (t) of 
the cavity output state \(fi). Suppose now that T is 
the final time of the interaction (i.e., the operation time 
determined by the driving laser pulse is from to T). 
The pulse shape / (t) is connected with the coefficients 
c u {£) before the frequency components in |</?i) by the 
Fourier transformation |29(| 



/(*) 



1 



'2-k J- 



dwc w (T)e- lw( *- T) . 



(17) 



From Eqs. (14), (16), and (17), we finally obtain 

/ (t) = sin 6 (t)exp(- J j sin 2 8 (t) dr^j . (18) 

Note that the single-photon pulse shape / (t) is com- 
pletely determined by 9(t), i.e., by the driving pulse 
shape e(t), and is independent of the random atom po- 
sition r and the absolute value of the coupling coefficient 
g(r). As we have mentioned before, this is the main 
advantage of this adiabatic method compared with ei- 
ther the Raman scheme or prior proposals based upon 
adiabatic passages with uniform illumination |25| , [l7| , 18| , 
and this feature is essential for many applications of this 
setup. 

The above result is obtained within the adiabatic ap- 
proximation, and in the adiabatic limit, the solution is 
independent of the atomic spontaneous emission rate 7 S 
and the detuning A. This is only a rough picture. In the 
following, we will solve exactly the dynamics governed 
by the Hamiltonian (6) without the use of the adiabatic 
approximation. The exact solution is necessary in the 
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following two senses: first, we need to verify the above 
ideal picture and to find out under what condition this 
picture is approximately valid. Though in the three-level 
case, we have some simple estimation of the condition 
for the adiabatic following, it is not easy to figure out 
the exact adiabatic following condition for more realistic 
situation of a continuum of external modes. In this case, 
the argument based on the level spacing is not valid. We 
need to know how long the operation time T should be to 
satisfy the adiabatic following condition. We also expect 
that the atomic spontaneous emission can not be made 
negligible simply by increasing the operation time T. Its 
rate j s should be small enough to satisfy the strong cou- 
pling condition Kjs -C q 2 , where ~g denotes the average of 
the coupling rate g (r) [p2| . Second, in real experiments, 
the operation time T is not infinitely long, and the cou- 
pling rate ~g can not be arbitrarily larger than the decay 
rates k and 7 S due to limitation of the technology (for in- 
stance, in Caltech experiments, typically, g/2n is around 
20MHz, and k/2tt ~ 7 s /2tt ~ 6MHz). In this case, there 
would be various non-adiabatic corrections to the above 
ideal picture, for instance, the atom may go down from 
the level |e) to |s) through a spontaneous emission, and 
then we lose the emitted photon and thus have no out- 
put from the cavity; or we have a single-photon output, 
but it is in a wrong and unknown pulse shape due to its 
sensitivity to the random atom position induced by the 
non-adiabatic contributions. It is desirous to calculate 
quantitatively the magnitudes of these noises to predict 
the real experiments. The exact solution of the system 
dynamics is only available with the numerical methods, 
which is the main task of the following section. 

III. EXACT NUMERICAL SIMULATIONS 
A. The numerical calculation method 

In this section, we solve exactly the system dynamics 
governed by the Hamiltonian (6) through numerical sim- 
ulations, and calculate various nonadiabatic corrections 
and noise magnitudes. For numerical simulations of the 
Hamiltonian (6), we need to discretize the free-space filed 
b (uj) by introducing a finite but small frequency interval 
Suj between two adjacent modes. Then, in total we have 
about N ss 2wb/5ui free-space modes, with the j mode 
denoted by bj . The frequency detuning ujj of the j mode 
is given by ujj = (j — N/2)5ui. To assure that there is 
no change of the physical result after the discretization, 
we should choose the frequency interval Suj much smaller 
than the inverse of the operation time T, and the band- 
width uib much larger the cavity decay rate k. 

For the numerical simulation, we can similarly expand 
the state \"^>) of the whole system in the form of Eq. (9), 
with the single-photon pulse state replaced by 



N 

IVi) | vac) . (19) 

3=1 

From the Hamiltonian (6), we get the following complete 
set of equations for the coefficients Cd, Cb, c e , and Cj 

N 

c d = - 6 c b — k' sin 8 ^ Cj , (20) 

3=1 

N 

c b =9 c d - i^/n 2 (r, t) + g 2 (r)c e + k' cos 8^ Cjl (21) 

3=1 



c e = HA - 7 ,/2) c e - i^/n 2 (r, t) + g 2 (r)c b , (22) 

Cj= —i (j — N/2) SuiCj + k' s'mOcd ~ k' cos 9cb, (23) 

where the effective decay rate k! = <J kSuj/2it. We ob- 
tain the solutions of the these coefficients by numerically 
integrating Eqs. (20)-(23) from the time t = to t = T, 
where T is the duration of the driving pulse e(t). We as- 
sume that e(t) is a Gaussian pulse so that a (t) is a Gaus- 
sian function of the time t, with its peak value at T/2, 
and a width t w significantly smaller than T/2. All the 
functions of 8 in Eqs. (20)-(23) are decided from cos9 = 

l/^Jl + \r a(t)\ 2 , and (r,i) + g 2 (r) = g (r) / cos 6. 
To simulate the randomness of the atom position r, we 
vary the value of g (r) in the simulation to look at whether 
the final result changes with this variation. 

B. Shape of the output single-photon pulse 

The output single-photon pulse shape / (t) can be eas- 
ily constructed from the solution of the coefficients Cj 
through a discrete version of Eq. (17). The result is 
shown in Fig. 2 for g (r) = 3k and g (r) = 6k. Although 
we have not made definitive measurements, wc estimate 
that g (r) varies within a factor of roughly 2 in the current 
Caltech experiment |^,^,^6l. Here and in the following, 
the pulse shape function fit) is always re-normalized ac- 
cording to J |/ (t)\ 2 dt = 1 for convenience of comparison. 
We see that the two curves overlap very well, which con- 
firms the prediction that the output pulse shape is very 
insensitive to the randomness of the coupling coefficient 
g (r) when the adiabatic condition is satisfied (we take 
T = 20/k for this figure). We also draw in this figure 
the pulse shape / (t) given by Eq. (18) derived in the 
ideal adiabatic limit, which agrees well with the exact 
numerical results. Therefore, within the adiabatic condi- 
tion, we can use the analytical result (18) to design the 
shape of the output single-photon pulse by modulating 
the driving pulse shape e(f). 
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FIG. 2. The shape of the output single-photon pulse de- 
scribed by the amplitude \f (t)\ versus the time t for the cou- 
pling rates g (r) = 3k (solid curve) and g (r) = 6k (dotted 
curve). The dashed curve represents the pulse shape in the 
ideal adiabatic limit calculated from Eq. (16). In this figure, 
we have taken 7 S = k, A = 0, and T = 20/k. The driving 
pulse e (t) is in a Gaussian shape with the peak at t = T/2 
and a width t w = T/5. 



C. Noise magnitudes and the adiabatic condition 

To quantify the noise magnitudes in this setup, we can 
define several error probabilities. First, we have the leak- 
age error due to the atomic spontaneous emission. A 
photon may be emitted to modes other than the princi- 
pal cavity mode through the spontaneous emission with 
the rate 7 S . As a result, the norm \a\ 2 + |c;,| 2 + |c e | 2 + 



we can use 



of the state (10) decays with the time t, and 



JY 



P spon = 1 - \ Cd (T)\ 2 - \c b (T)| 2 - |c e (T)| 2 - Y, \°i ( T )! 2 

3=1 

(24) 

at the final time T to quantify the total possibility of the 
spontaneous emission loss. Second, due to the finiteness 
of the operation time T and the pumping field amplitude 
e(t), the initial excitation in the dark state is not neces- 
sarily fully transferred to the output quantum signal at 
the final time, and we can use 



-PtranHcd (r)| 2 + | C 6 (T)f 



(25) 



at the time T to quantify the transmission inefficiency. 
In principle, we can arbitrarily decrease the transmission 
inefficiency by increasing the duration T or the amplitude 
e(t) of the pumping field. Finally, even if a photon is 
emitted into the cavity output field, it is not necessarily 
in the right pulse shape as given by Eq. (18) due to the 
non-adiabatic correction. This non-adiabatic correction 
depends on the random atom position and is unknown, 
so it is also a source of noise. To quantify this noise, we 



denote the ideal pulse shape given in Eq. (18) as fid (t), 
and the real pulse shape calculated from the numerical 
simulation as / rea i (t), then the shape mismatching error 
can be described by 



p . — 

± mis 



Jc r /real (*) /id (*)di 



Jo T |/rc al Ml^^l/id (t)\ 2 dt 



1/2 



(26) 



This quantity is directly related to the visibility of the 
fringes if we interfere two single-photon pulses from two 
such setups. 

For the example shown in Fig. 2, with g (r) = 3k = 37 s 
(the other parameters are given in the figure caption), we 
have P spon » 4.0%, P tran « 0.04%, P mis « 0.18%. The 
dominant source of noise is the leakage error P sp0 n in- 
duced by the spontaneous emission. If we increase the 
operation time T so that the adiabatic condition is bet- 
ter satisfied, the above-defined noise magnitudes can be 
reduced a little bit, but not too much. For instance, with 
the above example but T = 30/ k, we have P spon w 3.33% 
and P m i S w 0.15%. On the other hand, if T is reduced so 
that the adiabatic condition is not well satisfied, the error 
probabilities can significantly increase. Fig. 3 shows the 
output pulse shapes for g (r) = 3k and g (r) = 6k with 
T = 5/k. The two curves are obviously different to each 
other and are also different to the ideal shape as given 
by Eq. (18). For the example with g (r) = 3k = 37 s 
and T = 5/k, we have P sp0 n ~ 36%, P tr an ~ 3.2%, 
Pmis ~ 2.7%. All the noise magnitudes significantly 
increase. In particular, the spontaneous emission loss 
becomes very big. This can be easily understood since 
without the adiabatic condition, the excited state |e) will 
be populated during the operation, and thus we have a 
correspondingly larger spontaneous emission loss. 
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FIG. 3. The shape \f (t)\ of the output single- photon pulse 
for the coupling rates g (r) = 3k (solid curve) , g (r) = 6k 
(dotted curve), and in the ideal adiabatic limit (dashed curve). 
We assumed the same condition as in Fig. 2, except that 
T — 5/k, which does not well satisfy the adiabatic condition. 
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D. The strong coupling condition 

Next we look at the requirement of the strong coupling 
condition. Let g denote the average value of the coupling 
rate g (r) . Normally one requires ~g 2 ^> k^ s to satisfy 
the strong coupling condition. We can define the strong 
coupling parameter d sc as d sc = ~g 2 /Kj s , and calculate 
the above-defined noise magnitudes P S pon,Ptran,Pmis un- 
der different values of the parameter d sc . We assumed 
T = 30 /k and A = in the calculation so that the 
adiabatic condition is well satisfied. It turns out that 
the spontaneous emission loss P sp on is always the dom- 
inant loss (about ten times larger than other sources of 
noise). Thus, in Fig. 4, we only show the calculation 
result for P sp0 n under different values of d sc . The result 
can be approximately simulated by an empirical curve 
with P spon S3 1/ (4d gc ). 

We can use this simple formula to estimate the spon- 
taneous emission loss under different experimental con- 
ditions. Actually, in current experiments, the strong 
coupling condition is only marginally satisfied. For 
instance, for the cesium atom in the Caltech group, 
(k, 7s) /2tt ss (8,5.2) MHz (note that k and 7 S here de- 
note the energy decay rates, which are two times the 
corresponding amplitude decay rates) [|||], and g/2-K is 
expected to be approximately 15 MHz for the transition 



(6S 1/2 ,F = 4,m = +4) 



(6P 3/2 ,F = 4,m = +4) 



(Note that the transition (6S1/2, F — 4, m = +4) — > 
(6-P3/2, F — 5, m = +5) cannot be used as 
a A— configuration though it has a slightly larger cou- 
pling rate g). These values lead to d sc = g 2 /nj s S3 5.4 
and a resulting spontaneous emission loss around 4.6%, 
which is quite accessible with the present technology. As 
another example, in the recent experiment flSfl , one has 
(k,7 s ) /2tt w (1.25,6.0) MHz, and g/2ir S3 2.5 MHz ac- 
cording to the estimation there. With these parameters, 
we estimate that the spontaneous emission loss is about 
Pspon ~ 30% if one uses the scheme here. If the usual 
adiabatic scheme is adopted with a uniform driving pulse 
perpendicular to the cavity axis, the spontaneous emis- 
sion loss should be still significantly larger, as will be seen 
from the simulation in the last subsection. 



' spon 

0.15 



FIG. 4. 

average coupling rate g(in the unit of the cavity decay rate 
ft). We assumed A = and 7 S = «, so the strong coupling 
parameter d sc is simply (^/k) 2 . The circles represent the 
results from the numerical calculation, and the dashed curve 



is from the empirical formula P sl 
simulates the numerical results. 



1/ (4d sc ) which well 



E. The influence of the single-photon transition 
detuning 

In the above calculations, we assumed A = 0. Finally, 
we discuss the influence of a non-zero single-photon de- 
tuning A. In Fig. 5, we show the calculation result of 
the exact pulse shape function / rea i (t) with a significant 
detuning A = k, and compare it with the ideal pulse 
shape function /id (t) given by Eq. (18) for both of the 
amplitude and the phase. The other parameters for this 
example are given in the figure caption. From the figure, 
we see that the two amplitudes |/ rea i (t)\ and |/id (t) \ still 
overlap very well, but their phases become a bit different 
due to the detuning. 

This phase difference is determined by the the detun- 
ing A, whereas the latter depends on the different level 
shift between ground and excited states, and hence varies 
with atom position within the FORT beam. In the case 
of the simple level scheme depicted in Figure 1, the states 
\g) and |e) would have spatially dependent level shifts of 
opposite sign, which would lead to variations in A com- 
parable to the trap depth. Fortunately, there is a sim- 
ple avenue to mitigate this difficulty by considering the 
multi-levels involved for the FORT beam, as described 
in Ref . |33) , so that the trapping potentials for the states 
\g) and |e) are very nearly the same. For example, for 
the experiment of Ref. |lq| , the difference in trap depth 
for \g) and |e) is roughly 10% of the trap depth. Relative 
to the current analysis, there is then a variation in A 
as the atom moves in the FORT potential, which is un- 
known when the adiabatic protocol is implemented. The 
curve in Fig. 5 is an attempt to estimate the impact of 
such random detunings by setting A = k, which exceeds 
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the actual magnitude of any spatially dependent detun- 
ings for FORT depths up to about 50MHz. The phase 
difference in the pulse shape function caused by the un- 
known detunings is a source of noise, which contributes 
to the shape mismatching error defined in Eq. (26). For 
this example with g (r) = 3k, we have P sp0 n ~ 3.33%, 
-Ptran ~ 1CP 4 , which are basically the same to the cor- 
responding case without detuning, but P m i s w 3.33%, 
which becomes significantly larger due to the contribu- 
tion of the phase difference. 

0.2 1 1 1 1 1 1 1 




5 10 15 20 25 30 

t 

FIG. 5. The amplitude (the modular) and the phase (di- 
vided by 7r/2) of the real pulse shape / rea i (t) (two solid curves) 
and the ideal pulse shape /id (t) (two dashed curves) versus 
the time t with the single-photon transition detuning A — k. 
We assumed g (r) = 3k, *y 3 — k, and T = 30/k. In this 
case, the main difference between / rea i (t) and /id (t) lies in 
the phase difference. 

F. Compare with the usual adiabatic scheme 

In our scheme, the driving pulse is matched to a cavity 
mode which has basically the same spatial mode struc- 
ture as the cavity QED light. In usual adiabatic schemes 
p5j^8| , the driving laser is assumed to be perpendicu- 
lar to the cavity axis with uniform illumination intensity. 
We expect that with the present interaction configura- 
tion, our scheme is more insensitive to the randomness 
in the atom position. To compare the two configurations 
more quantitatively, we have calculated the output pulse 
shapes and noise magnitudes for both schemes. 

First, let us assume that the atom has been trapped in 
one potential well, but the coupling rate g (r) may vary 
within a factor 2 due to the unknown atom position. 
In figure 6, we show the calculation results of the out- 
put pulse shapes. The solid curve shows the pulse shape 
function \f (t)\ when g (r) = 3k and fi TO (r) = 3k, where 
Q m (r) is the maximum of f2 (r, t) with respect to time t 
(f2 (r, t) is assumed to be a Gaussian function of i as spec- 
ified in the caption of Fig. 2). Now, if g (r) varies by a 
factor 2 due to change of the atom position, in our scheme 



the Rabi frequency will correspondingly change by the 
same ratio. The dashed curve shows the pulse shape 
for g (r) = 6k and Q m (r) = 6k. One can see that the 
two curves overlap very well with the mode mismatching 
noise smaller than 0.2%. In contrast, in usual adiabatic 
schemes with uniform illumination intensity, f2 m (r) does 
not change as g (r) varies with the atom position, so we 
have the same f2 m (r) = 3k. The dotted curve in Fig. 6 
shows the pulse shape for g (r) = 6k and f2 m (r) = 3k. 
It is significantly different to the above two curves with 
a notable mode mismatching noise P m i s ss 6.9%. There- 
fore, by this interaction configuration, the scheme is more 
robust to the random variation of the atom position. 
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FIG. 6. The shape |/ (t)| of the output single- photon pulse 
for the following pairs of coupling rates and the driving Rabi 
frequencies: first, g (r) = 3k and Q m (r) = 3k (solid curve), 
second, g (r) = 6k and fi m (r) = 6K(dashed curve), and fi- 
nally, g (r) = 6k and O m (r) = 3k (dotted curve). The other 
parameters are the same as in Fig. 2. 

The improvement by this protocol is more remarkable 
if we consider the case where the atom is not fixed in 
one potential well, and may jump form well to well in 
the axial direction, as is the case in some recent exper- 
iments fl^ , |l8| . The variation of the atom position in 
the axial direction is typically fast compared with the 
operation time T, so we have a time-varying atom po- 
sition r and coupling rate g (r) . Here, we consider an 
explicit form of the time variation of g (r) by assuming 
g (r) = 6k sin (Ant/T + ipo), where the phase (po is ran- 
domly chosen corresponding to the randomness in the 
initial atom position. It is enough to illustrate the gen- 
eral result by considering this special example. First, let 
us calculate the output pulse shape / (t) for the usual 
adiabatic scheme where TO (r) is fixed as a constant 
p5|,^7|Jl8|. The solid and the dashdot curves in Fig. 7 
show the real parts of / (t) with initial phase tpo = and 
tpo = 7r/2, respectively (the imaginary parts of / (t) are 
actually small and negligible). The two curves do not 
overlap at all. Neither the magnitude nor the phase of 
the pulse shape / (t) can be controlled with this scheme. 
We also calculate the spontaneous emission loss P sp0 n for 
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this example. The average spontaneous emission loss is 
about P sp on ~ 25%. 
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FIG. 7. The real part of the pulse shape func- 
tion Re [f (t)] as g (r) varies with time by the form 
g (r) = 6k sin (4nt/T + ipo) in the usual adiabatic scheme 
with ifo = (solid curve) and </?o = 7r/2 (dashdot curve), 
respectively. The dashed curve shows the ideal pulse shape 
calculated from Eq. (18). The other parameters in this figure 
are the same as those in Fig. 2. 

Similarly, we can calculate the pulse shape for the same 
example with the present scheme. In this case, due to the 
atomic motion, fl m (r) varies with time in the same way 
as g (r), but the ratio f2 m (r) j g (r) is kept constant. Fig. 
8 shows the real part of the shape function / (t) in this 
case, with the solid and the dashdot curves corresponding 
to the initial phase ipo = and ipo — 7r/2, respectively. 
Although the two curves do not overlap very well, they 
still look similar with the same phase. They also roughly 
agree with the ideal shape function given by Eq. (18), 
which is shown as the dashed curve in Fig. 8. The aver- 
age mode mismatching noise for these two curves is given 
by -Pmis ~ 1.1%, and the average spontaneous emission 
loss is P sp0 n ~ 9.4%. The spontaneous emission loss is 
also significantly reduced with the present scheme. This 
can be understood as follows: if one has a constant O m (r) 
as the usual adiabatic scheme, when the atom moves to 
the place with g (r) near to zero, the adiabatic condition 
is not well satisfied, and as a result, one has a consid- 
erably large spontaneous emission loss; however, in the 
present scheme, in the place where g (r) is near to zero, 
Q m (r) is also near to zero. The excitation probability 
of the atom is then reduced, and the adiabatic condi- 
tion is better satisfied. Consequently, one has a smaller 
spontaneous emission loss. 
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FIG. 8. The real part of the pulse shape function Re [f (t)] 
calculated for the same example as in Fig.7, but now for the 
present adiabatic scheme where the driving Rabi frequency 
varies in the same way as the the coupling rate when the 
atom moves. 



IV. APPLICATIONS 

There have been many proposals to use the setup 
with single atoms in high-Q cavities for various appli- 
cations, such as for single-photon or entangled photon 
source |(|f7j , for quantum communication between differ- 
ent cavities || , for atomic quantum teleportation jll],[l2| . 
and for quantum computation jl3| . In these proposals, 
one always assumed the atom is well localized so that 
the Lamb-Dick condition is satisfied. However, one can 
apply the method here to all of the schemes mentioned 
above to eliminate the challenging Lamb-Dick condition. 
Basically, what one needs to do is to replace the Raman 
scheme with the adiabatic scheme, and to keep the pump- 
ing laser collinear with the cavity axis so that the driving 
pulse and the cavity mode have the same spatial mode 
structure. All the calculation results (for the noise mag- 
nitudes, pulse shape, etc.) in this paper apply to these 
schemes. After the improvement, it becomes consider- 
ably easier to implement these schemes with the current 
technology. Here, we briefly mention some possibilities. 

A. Controllable single-photon or entangled-photon 
source 

It is desirable to have a single-photon source with 
all its properties fully controllable, including its emis- 
sion direction, emission time, and pulse shape. This 
kind of source has important application in some recent 
quantum information processing schemes [Q , which are 
normally based on the interference of different single- 
photon pulses. To get interference between deferent 
single-photon pulses, one needs to require all the pulses to 
be directional and to have the same time shape. Recently, 
there have been significant experimental advances in the 
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realization of the single-photon source |35| |3J7].|_L7j. 18 1 . In 

the experiments based on the solid-state material [ 35 -[37) , 
the single-photon emitter has a fixed position, and one 
can in principle well control the pulse shape. However, 
the emitted pulse is typically not directional. On the 
other hand, in current experiments 0Jl| with high-Q 
cavities, the emitted pulse is directional, but its shape is 
not well controlled since with uniform illumination of a 
perpendicular driving pulse, the waveform f(t) depends 
on the time history of the coupling rate g (r), which in 
turn depends on the atom's position. As the atom falls 
through the cavity, it has basically a random trajectory, 
leading to unknown variations in g (r) both in magni- 
tude and sign. It is a challenging experimental endeavor 
to demonstrate a single-photon source with all the prop- 
erties mentioned above fully controllable. 

The method in this paper shows that the single atom 
trapped in a high-Q cavity is a good candidate for the 
realization of the fully controllable single-photon source. 
Though the coupling rate g (r) is not completely fixed in 
current setups due to the hardness to fully localize the 
atom, the emitted single-photon pulse has a definitely 
well controllable time shape and emission direction with 
an appropriate design of the interaction configuration as 
has been shown before. 

As shown in Ref . , with a more involved atomic level 
structure, it is possible to engineer entanglement between 
different single-photon pulses. It is straightforward to 
combine the method here with that scheme to eliminate 
the requirement of the Lamb-Dick condition in so that 
one can get an entangled single-photon source with the 
"hot" trapped atom as well. 



B. Quantum communication between different 
cavities 

The dynamics governed by the Hamiltonian (6) is re- 
versible if we neglect the atomic spontaneous emission 
7 S . Therefore, if one directs the emitted single-photon 
pulse back to the cavity, and at the same time reverse 
both of the time shapes of the single-photon pulse and 
the driving pulse, the single-photon pulse will be com- 
pletely absorbed as long as the noise effects are negligi- 
ble. It was first proposed in || that one can use this 
kind of phenomenon to achieve quantum communication 
between different cavities, that is, to transfer quantum 
states of a trapped atom from one cavity to another cav- 
ity. For this purpose, one can require that the emitted 
single-photon pulse has a time-symmetric shape by mod- 
ulating the driving pulse shape. For a time-symmetric 
pulse, its time reversal is itself, so we can directly direct 
this pulse to another cavity with the same configuration 
but with a time-reversed driving pulse, then the single- 
photon pulse will be completely absorbed by this cavity, 
which transfers the atomic state from one cavity to the 
other one. The scheme in || is based on the Raman con- 
figuration, but it is straightforward to transfer it to the 



adiabatic configuration discussed in this paper so that 
it works with a hot trapped atom. Note that the same 
setup can also be used for storage of a single-photon pulse 
with a known shape j2^]3^|3^] . 

To get a time-symmetric single-photon pulse for a com- 
plete absorption of the second cavity, Ref. || gives a nu- 
merical solution to the shape of the driving pulse. For 
the adiabatic configuration, within a good approxima- 
tion, we have an analytic expression (18) which connects 
the shape of the output single-photon pulse to the shape 
of the driving pulse. With this analytical expression, it 
becomes easier to design the shape e(t) of the driving 
pulse. For this purpose, we solve sin#(i) from Eq. (18) 
as 



^sin^t) 



/(*) 



l-/J/ a (r)dr 



(27) 



The form of sin 9 (t) is immediately available from this 
equation for any desirable output pulse shape / (t) (which 
has been assumed to be real and positive for simplicity) . 
Then, the shape of the driving pulse can be easily decided 

from e(t) cc a (t) and sin 9 (t) = r Q a (t) / \Jl+ \r a a (t)\ 2 , 
where r a is the ratio of the Clebsch-Gordan coeffi- 
cients. For instance, if we want to have a time sym- 
metric / (t) in the period < t < T with the form 
/(*) = y/p/2sech[P(t~T/2)} , where we have as- 
sumed sech (— (3T/2) <C 1, sin9 (t) should be in the form 

»(*) 



sin V (t) = yJP/KsJl + tanh \fi (t - T/2)}. Note that we 
only have a solution of 9 (t) when the rate < k/2, 
which is consistent with the observation that any pulse 
from the decay of a cavity can not vary with the time 
faster than the cavity decay rate. From sin0(i), we see 
that the shape e(t) of the driving pulse should be chosen 
according to 



e (t) oc 



1 + tanh [/3 (t - T/2)] 
{k/(3 - 1) - tanh \p (t-T/2)]' 



(28) 



As a special case, if k//3 = 2, e(t) oc e l3t ^ t ~ T / 2 \ which 
grows exponentially with the time t for the operation pe- 
riod < t < T. Therefore, we have a simple solution to 
the driving pulse shape for quantum communication be- 
tween two different cavities: for the first cavity, we apply 
an exponentially increasing pulse with e(t) = e(0) e Kt / 2 , 
and for the second cavity we apply its time reversal, 
that is, an exponentially decreasing pulse with the decay 
rate k/2. The pulse duration T should satisfy kT ^> 1, 
and the initial value e (0) is determined by the require- 
ment r a a {T/2) = 1. The single-photon pulse connecting 
the two cavities then has a time symmetric shape with 
f(t) oc sech [K(t- T/2) /2]. 



11 



C. Entanglement generation and atomic quantum 
teleportation 

If one has two cavities, each with an atom inside, one 
can maximally entangle these two atoms 1 and 2 by the 
following method: The two atoms are initially prepared 
in the state \g), and then we excite them to the state \s) 

with a small possibility po sa 1— exp Jq sin 2 (r) dr^j 

through an incomplete adiabatic passage. The output 
pulses from the two cavities, each with a mean photon 
number po , have a definite pulse shape as we have shown 
before so that they can interfere with each other at a 
50% — 50% beam splitter. The outputs of the beam 
splitter are detected by two single-photon detectors, and 
if we register a photon from one of the detectors, due 
to the interference, we do not know from which cavity 
the registered photon comes from. The two atoms 1 and 
2 are thus projected to a quantum superposition state 
(15)1 l s )2 ^ l s )i 15)2) /V2) which is maximally entangled. 
The method described here is just an adiabatic passage 
version of the scheme in Refs. pT,^2|. By transformation 
from the Raman version to the adiabatic passage version, 
the output pulse shapes become insensitive to the ran- 
dom atom position as is required for interference, which 
is important for the scheme to work with hot atoms. 

After entanglement has been generated, one can use 
it for atomic Bell inequality detection, for quantum tele- 
portation of atomic states or even for realization of 
quantum repeaters jf0). To realize quantum repeaters, 
what one needs to do is to simply replace the atomic en- 
semble in the scheme in Ref. |4C| by the setup of a single 
atom in a high-Q cavity. 

For the above applications, in addition to the entan- 
glement generation, we also need to do some single-bit 
operations. These single-bit operations should also be 
performed in a suitable way so that they are insensitive 
to the random atom position r. One way is to still use 
adiabatic passages. It is possible to realize any single-bit 
operation with adiabatic passages but for this 

purpose one needs to use a four-level scheme instead of 
the A configuration. There is actually a simpler way 
for getting robust single-bit operations based on the Ra- 
man transitions. Note that for single-bit operations, we 
do not need to use any cavity mode or cavity effect. 
We can shine two travelling-wave beams on the atom 
respectively coupling to the transitions \g) — > |e) and 
\s) — > |e). They are assumed to be collinear and prop- 
agating along the x axis which is perpendicular to the 
cavity axis z. The two travelling- wave beams are broad 
with the beam radius much lager than the typical vari- 
ation length of the atom position. With this condition, 
the two Rabi frequencies for the transitions \g) — > |e) 
and \s) — > |e) are given by Qi (r) = flioe lu>9eX ^ c and 
f) 2 (r) = f^oe 1 "" 2 ^, respectively, where fiio and f2 2 o 
are basically independent of the atom position r. Under 
a lager detuning A, the effective Raman coupling rate 
flu. ~ fl\ (r) £1 2 (r) / A oc e m,,I ' c , which is very insensi- 



tive to the random atom position r since cjuj gs is typically 
much larger the variation length of the position. There- 
fore, as long as we do not need to use the cavity effect, 
a Raman scheme with two broad collinearly-propagating 
beams suffice to eliminate the sensitivity to the random 
atom position. 



D. Quantum computation 

In principle, we can also use this setup for quan- 
tum computation jlq] , and eliminate the requirement of 
the Lamb-Dick condition by performing all the quantum 
gates using adiabatic passages |4l],^| with appropriate 
configurations. However, the requirements for a univer- 
sal quantum computation are more challenging compared 
with the applications mentioned above, and this it is 
somewhat a long-term goal, so we do not discuss here 
the details of this possibility. 



V. DISCUSSION OF THE EXPERIMENTAL 
SITUATION 

Finally, let us mention the current experimental situa- 
tion related to this work at the Caltech group. In the Cal- 
tech experiment, a single cesium atom is trapped inside 
the high-finesse cavity with a FORT beam. The atomic 
states |<?), |s), and |e) correspond to the hyperfine lev- 
els (65i /2 , F = 3, m = +3) , (651/2, F = 4, m = +4) , and 
(6P3/2, F = 4, m = +4), respectively. The FORT beam 
is incident on one of the cavity mirrors and resonant to 
a longitudinal mode of the cavity. Presently, the FORT 
wavelength Afort is 936 nm. This wavelength was cho- 
sen because with such a beam, the trapping potentials for 
the ground 6<S f 1 / 2 manifold and the excited 6P 3 / 2 mani- 
fold are nearly identical. Considering only this reduced 
manifold of states, we find that the expression for the 
FORT potential of the ground states \g) and \s) is given 
by @ 

Ufort(v) = ^-?(-r- + -r-)I(r). (29) 
2luq A 2 Ai 

Here, Ai(A 2 ) is the detuning of the FORT light of fre- 
quency wfort = 27rc/ Afort from the Pi/ 2 (P 3/2 ) level, 
and 7 s /27r w 5.2 MHz is the spontaneous decay rate of 
the level 6P 3 / 2 . The intensity J(r) of the standing wave 
mode inside the cavity is given by 



J(r) 



8P 



■sin 2 ( 



2ttz 
Afort 



) exp(- 



(30) 



where Wg ~ 25 /im is the waist of the Gaussian mode, 
and P is the power of the FORT beam inside the cavity. 
The trap frequencies faxiai, ^radial in the axial and radial 
directions follow from these expressions as 
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(*Wl^adxal) = J 2U ^OBT 1 2 rj -^-), 

Zirn y mc y m{wo) 

(31) 

where Uq = Ufort{®) is the trap depth. The typical 
power of the FORT beam measured outside the cavity is 
about 1 mW, and the power P inside the cavity is en- 
hanced by a factor of the cavity finesse, which is about 
2200 at the wavelength of the FORT beam. With this 
number, the typical values for the trap depth and fre- 
quencies are given by Uq ss 38 MHz, ^ a xiai ~ 510 kHz, 
and Radial ~ 4.3 kHz, respectively. The current achiev- 
able temperature T tcm of the trapped atom is a significant 
fraction of the trap depth Uq (such as a half). With such 
a temperature, the spatial extent of the atomic motion in 
the axial and radial directions are estimated respectively 
by 

Sz/ Afort ~ (1/271") arcsin yj k B T tcm /U , (32) 



6rj_ w w y/- ln(l - k B T tcm /U ), (33) 

which will induce significant variation of the coupling rate 
g (r) given by Eq. (1). For example, for the temperature 
of half of the trap depth, the axial uncertainty is 120 nm, 
while the radial one is 15 /im. These uncertainties cause 
variations in g of 30% due to the radial motion, and 35% 
due to the axial one. Therefore, within the current ex- 
perimental technique, it is important to use the method 
given in this paper to make the application schemes in- 
sensitive to the variation of g (r) . The time scale for 
the variation of g (r) is estimated by the inverse of the 
trap frequencies z^axiai, ^radial in the axial and radial di- 
rections, respectively. The operation time T is typically 
significantly shorter than 1/V r adiai> but longer or compa- 
rable to 1/^axiai- So, we can take static average of g (r) 
in the radial direction, and dynamical average of g (r) in 
the axial direction as discussed in the note ||32| . 

We also would like to note that although the method 
in this paper shows that many application schemes of 
the cavity QED setup can be demonstrated before the 
achievement of efficient cooling of the trapped atom in- 
side the cavity, the cooling is still an important and de- 
sirable technology yet to achieve to significantly increase 
the trapping time of the atom. In addition, a combina- 
tion of the cooling technology and the method here could 
further improve the performance of various application 
schemes. 



VI. SUMMARY 

In summary, we have shown that the setup with a sin- 
gle trapped atom in a high-Q cavity can be used to re- 
alize many cavity QED and quantum information pro- 
cessing schemes even if the atom is still hot and not fully 



localized in space (the Lamb-Dick condition is not yet 
satisfied). This could significantly simplify the on-going 
experiments since it means many interesting schemes can 
be demonstrated with the present technology before the 
achievement of efficient cooling inside the cavity. Even 
with further advances in atomic localization in cavity 
QED, our scheme should lead to a greater robustness 
again certain experimental nonidealities. The basic idea 
of this method is to design an appropriate adiabatic pas- 
sage so that the relevant dynamics only depend on the 
ratio of two coupling rates. Though each of the cou- 
pling rates is sensitive to the unknown or time-varying 
atom position, their ratio is fixed and controllable as 
the two rates depend on the random atom position in 
the same way with the appropriate interaction configu- 
ration that we have described. We confirm the valid- 
ity of this method by solving the complete model which 
describes the realistic setup. The approximate analyt- 
ical solution and the exact numerical simulations agree 
with each other. From the numerical simulations, we also 
calculate quantitatively various noise magnitudes in this 
setup, and show one can achieve reasonably good perfor- 
mance with the values of the parameters based on the 
present technology. Finally, we show that this method 
can be incorporated into many previous schemes, allow- 
ing the demonstration of these application schemes with- 
out the requirement of the full localization of the atom. 
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